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We propose a new approach to study the evolution of a quantum state that is 
encoded in a system which is continuously subject to the operations required 
to implement a quantum error correcting code. In the limit of continuous error 
Q\ ', correction we introduce a Markovian master equation that includes the effects of: 

a) Hamiltonian evolution, b) errors caused by the interaction with an environment 
and c) error-correcting operations. The master equation is formally presented for 
all stabilizer codes and its solution is analyzed for the simplest such code. 

o 

On 
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1. Introduction 

-i— > ' 

In the last two years there has been enormous progress in the development of 
quantum error correction codes. The basic goal of these techniques is to protect 
quantum information against the degrading effect of uncontrolled interactions 
with external agents. To achieve this goal, being able to partially restore the 
unknown quantum state, one should perform a sequence of operations on the 
system used to represent the information. This active process, a quantum error 
correcting code, involves the following basic ingredients: First, one should encode 
the information (initially contained in the quantum state of k qubits) into an 
entangled quantum state of n carriers (n > k). After the encoding, the original 
state is protected from certain errors provided one acts on the system applying 
a so-called recovery operator. The simplest (but not the only) way to describe 
this operation is to decompose it in the following steps: one should first perform 
a measurement on n — k qubits, whose result reveals a syndrome. This measure- 
ment is followed by a unitary operation on the remaining k qubits, whose nature 
depends upon the syndrome. Finally, the operation is completed by reseting the 
n — k measured qubits to a reference (ground) state. The short history of quantum 
error correcting techniques is quite remarkable: after the first pioneering papers 
by P. Shor and A. Steane (Shor, 1995; Steane 1996a, b), several papers helped 
in clarifying the nature of quantum codes (Laflamme et at, 1996; Bennett et 
al, 1996). Very soon, the mathematical structure behind these simple physical 
ideas was revealed (Knill & Laflamme, 1997) and the, rather general, theory of 
stabilizer codes was fully developed (Gottesman, 1996a, 1997; Calderbank et al, 
1996). More recently, it was shown that the use of concatenated quantum error 
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correcting codes (Knill & Laflamme, 1996) together with fault tolerant techniques 
for quantum computation (Shor, 1997) may allow arbitrarily long computation 
provided the accuracy per operation is below certain computational threshold 
(Knill, Laflamme & Zurek, 1997; Preskill, 1997). 

In this paper we propose a new tool that could be used to study the evolution 
of an encoded quantum state. The motivation of our method lies on the inter- 
esting analogy between the active mechanism of quantum error correction and 
the process of cooling a system pumping energy away from it by driving certain 
internal transitions. Thus, we propose to use a master equation (inspired in the 
quantum-optical master equation) to describe the limit in which a system formed 
by n-carrier qubits is continuously being affected by the active process of error 
correction while, at the same time, is suffering the continuous interaction with an 
environment and, eventually, is evolving under a specified Hamiltonian. The limit 
of continuous error correction is the one in which one assumes that all operation 
involved in the code (encoding, applying the recovery operator, refreshing the 
syndrome, etc) can be considered to be instantaneous for all practical purposes. 

Our purpose here is twofold. On the one hand we will introduce the general 
master equation for a rather large class of error correcting codes (stabilizer codes) . 
On the other hand, as we will show that this master equation can be explicitly 
solved in some simple but relevant cases, we will discuss important features of the 
solution. We find it rather interesting that quantum error correcting codes admit 
a description in terms of a continuous differential equation. But it is certainly 
most remarkable that this equation not only can be formally written in terms 
of some abstract operators but can also be explicitly solved in some interesting 
cases. The rest of the paper is organized as follows: In Section 2 we present a 
very brief introduction to the theory of stabilizer error correcting codes. The 
ideas discussed in this section enable us to introduce the notion of "corrective 
quantum jumps". In Section 3, this notion is used to motivate a master equation 
describing the evolution of a quantum system which is continuously being affected 
by errors and is subject to a continuous sequence of corrective operations. We 
analyze and solve the master equation for the simplest example of a stabilizer 
code: the n = 3 qubit code protecting k = 1 qubit of information against phase 
errors affecting any one of the three carriers. In this case we present the explicit 
solution of the master equation computing the fidelity from the density matrix 
of the system. We discuss the various time scales that affect the evolution of the 
encoded state in the collective Bloch sphere. 



2. Stabilizer error correcting codes. 

We will consider a quantum error correcting code that protects k qubits by 
encoding them into n carriers. The code space TL^ is a 2 k dimensional subspace 
of the Hilbert space of the n carriers. H n is a tensor product of n 2-dimensional 
factors and has a natural basis whose elements are product states of the individ- 
ual carriers. This "physical basis" can be formed with the common eigenstates of 
the operators {Z±, . . . , Z n } pQ, Yi and Zi denote the Pauli matrices for the i-th 
qubit). We will label states of this basis not by the eigenvalues of the correspond- 
ing operators (which are ±1) but by the eigenvalues of the projectors onto the — 1 
subspace (which are or 1): Thus, the label Zj = (zj = 1) corresponds to a +1 
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(— 1) eigenvalue of the operator Zj. Furthermore, we order the n carriers in such 
a way that the last k qubits are the ones whose state we encode while the first 
n — k are the ancillary carriers. Therefore, states of the physical basis are of the 
form \s, z) p = \s) p (g> \z) p (where the strings s = (s\, . . . , s ra _fc), z = (z±, . . . , z^) 
store the corresponding eigenvalues and the subscript P is used to identify states 
of the physical basis). 

An error correcting code is a mapping from the physical product states |0)p <g> 
\^) p onto the code space 7ik, which is formed by entangled states of the n carriers. 
A rather general class of codes can be described in terms of their stabilizer group 
(see Gottesman, 1997). The stabilizer of the code is an Abelian group formed 
by all operators which are tensor products of Pauli matrices and have TLk as 
an eigenspace with eigenvalue +1. Every element of the stabilizer, which is a 
finite group with 2 n ~ k elements, can be obtained by appropriately multiplying 
n — k generators, which will be denoted as M±, . . . , M n _fc. The elements of the 
stabilizer are completely degenerate in the code space Ti^ (since all states in Ti^ 
are eigenstates with eigenvalue +1 of all Mj). To define a basis in the code space 
we choose k operators Li, . . . , Lj. which being tensor products of Pauli matrices 
commute with all elements of the stabilizer. These operators Lj>, f = l,...,k 
are the "logical pointers" since they define the directions in Tik associated with 
the logical states |0 >l, . . . , |2 fc — 1 >l (logical pointers belong to the group of 
operators which commute with the stabilizer, known as the normalizer). 

The n — k generators of the stabilizer together with the k logical pointers are 
a Complete Set of Commuting Operators (CSCO) whose common eigenstates 
form a complete basis of the Hilbert space 7i n . Elements of this "logical basis", 
labeled by n quantum numbers, are denoted as \m, where the bit strings 
m = (mi, . . . , m n _fc), / = (l±, . . . , 1^) identify the corresponding eigenvalues and 
the subscript L refers to logical states. The CSCO formed by the generators of 
the stabilizer and the logical pointers defines a prescription for decomposing the 
original Hilbert space of the n-carriers into a tensor product of a 2 fc -dimensional 
logical space A and a 2 ra_fc -dimensional syndrome space S. In fact, elements of the 
logical basis (which are entangled states of the n-carriers) are tensor products of 
states belonging to A and S: |m, 1)l = |m) l ® Encoded states, which belong 
to Hk, are also product states of the form |^) = |0)l <8> J2i c i\1)l- 

The code protects quantum states against any error E a whose action on states 
of the logical basis is to change the logical syndrome and, eventually, rotate the 
logical state in A in a syndrome dependent way: 

E a \m) L ® \l) L = e*™ |m + c a ) L U a \l) L . (2.1) 

Here U a is a unitary operator acting on the collective logical space A and <p ma is a 
phase that may depend on the syndrome and the error. The error E a changes the 
syndrome from m to m + c a where c a is the bit string storing the commutation 
pattern between the error and the generators of the stabilizer (the j-th bit of 
this string is one if the error anti-commutes with Mj and is zero otherwise). The 
reason for this is that when acting on a logical state, the error E a changes the 
eigenvalue of the operator Mj only if {Mj,E a } = 0. The label a used to identify 
errors is arbitrary and, for the case of non-degenerate codes (which is the only 
one we will consider here) it is always possible to label errors E a using simply 
the commutation pattern c a (i.e., we can choose a = c a ). 

To correct against the action of any of the errors E a (or against any linear 
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superposition of them) one can first detect the error by measuring the collec- 
tive syndrome (i.e., measuring the observables Mj, j = 1, . . . , n — k) and later 
recover from the error by applying the corresponding operator U\ . This detection- 
recovery process can be conveniently described as a quantum operation defined 
by the following mapping from the erroneous density matrix pi n into the corrected 
one p out : 

N 

Pout = X! RmPinR m (2-2) 
m=0 

where the sum runs over all syndromes (N = 2 n ~ k — 1) and the recovery operator 
is 

R m = \0)l lH (8>C4. (2.3) 

By construction, these operators satisfy the identity J2m=o ^m^m = I- 

As our description of error detection-recovery process is entirely formulated in 
the logical basis it does not involve any reference to encoding-decoding operations 
which can be simply defined as a change of basis: The encoding operator C 
is a unitary operator mapping the physical basis, formed by product states of 
the n carriers, onto the logical basis, formed by entangled states. Accordingly, 
C transforms the operators Zi (whose eigenvalues define states in the physical 
basis) into the operators Mj, Ly (that label states in the logical basis). Thus, the 
encoding operator C is such that Zj = C^MjC, j = 1, . . . , n — k and Z n _k+j> = 
C^Lj/C, j' = 1, . . . , k. Taking this into account the action of the operator R m 
can be described, in the physical basis, as the following sequence of operations: 
i) decode the state, ii) measure the syndrome in the physical basis by measuring 
Zj in the first (n — k) carriers, iii) If the result of the measurement is the string 
s, apply the syndrome dependent recovery operator UJ resetting the syndrome 
back to zero, iv) encode the resulting state. 

Finding a stabilizer code correcting a given set of errors is a rather hard task 
which involves designing generators having appropriate commutation patterns 
with the errors. Once the generators are found and the logical pointers are cho- 
sen an encoding-decoding operator can be constructed (strategies for designing 
encoding-decoding circuits from the stabilizer are known; see Gottesman, 1997). 
The recovery operators depend on the encoding-decoding strategy and can be 
explicitly found from the encoding circuit by running errors through it. 

In the next Section we will illustrate our model for continuous error correction 
using the simplest stabilizer code: the one encoding k = 1 qubit using n = 3 
carriers correcting against phase errors in any of the carriers. The basic errors 
are E\ = Z±, E<i = Z2 and £3 = Z3. The stabilizer of the code can be chosen 
to be generated by: Mi = X1X3 and M2 = X2X3. The commutation pattern 
associated with each error is c\ = 01 (because the error Z\ commutes with M\ 
and anti-commutes with M2), C2 = 10, C3 = 11 (note that we are labeling errors 
with the commutation pattern). An encoding circuit is shown in Figure 1. The 
operator C has the following simple properties: 

C*Z 1 C = X 1 , C^Z 2 C = X 2 , C ] Z 3 C = X X X 2 Z Z . (2.4) 

These properties entirely determine the action of the errors Zi in the logical basis. 
For example, the last identity implies that E^\m)L\l)L = (—^) 1 \™+c?i)l\1)l- Thus, 
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Figure 1. Encoding circuit for the three qubit code protecting against phase errors in the physical 
basis. Qubits are labeled from bottom to top. The third qubit is the one to be encoded. When run 
from left to right the circuit defines the decoding operator D. The encoding operator corresponds 
to the same circuit run backwards. 

the error E% not only changes the syndrome but also modifies the logical state by 
adding a phase. This means that the recovery operator for this error is U3 = Z 
Analogously, we can find how the other errors act on the logical basis showing 
that Ui = U 2 = L 



3. Master equation with corrective quantum jumps. 

Quantum error correction is an active process: for the density matrix to evolve 
according to the rule ( |2.2[ ), a sequence of physical operations must be carried 
out on the system (measuring the syndrome, refreshing it and applying the syn- 
drome dependent recovery operator). Let us suppose that we are able to perform 
detection-recovery operations arbitrarily fast (i.e., the time required for these 
operations is negligible with respect to any other time scale in the system) and 
that the corrective protocol is repeatedly applied many times. In the limit of 
continuous error correction we will consider that the time required to perform 
the corrective operations goes to zero while the number of corrective events goes 
to infinity. In this limit the evolution of the density matrix is generated by the 



continuous application of the mapping (|2.2|) . 



In this paper we suggest that this rather idealized scenario, which might be 
useful to analyze the optimal protection achievable with a given code, could be 
described using a simple master equation. To motivate the master equation one 
could noti ce th at the mapping from the erroneous to the corrected density matrix 
(given by ( |2.2[) ) is naturally described in the language of quantum jumps: In fact, 
the operators R m define "corrective jumps" which reset the logical syndrome to 
m = while acting as in the logical space A. Thus, to describe the continuous 
limit we propose to use a Markovian master equation (analogous to the quantum 
optical master equation, see for example Walls & Milburn, 1994) having N + 1 
decay channels associated with the i? m -quantum jump operators: 

p = L 0P + Ira (R m pR ] m ~ -R^RmP ~ ^pR^R^ ■ (3.1) 

Here, the constants j m measure the strength of the corrective jumps and Lqp 
includes all other sources of evolution (produced by a Hamiltonian Hq or caused 
by the errors E a ). To write this non-corrective piece of the master equation we will 
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consider a system with a Hamiltonian Hq that is interacting with an environment. 
We assume that this interaction produces a term in the master equation which, 
in the Markovian limit, has the same form than (|3.1|) with N decay channels E a : 



N / i i \ 

L oP = -i[H , p] + Y,<{ E apEl - ^E\E aP - -pE\E a \ . (3.2) 



0=1 



The complete master equation can be written in a rather compact form if we 
assume that all corrective jumps have equal strength ( r y m = 7) and that the same 
is true for the erroneous jumps (7^ = 7'). In this case, and if the erroneous jump 
operator are unitary, the complete master equation reads: 



N N 



p = -%[H , p) - (7 + Ni)p + 7 £ RrnpRl + l' £ E aP E l ( 3 ' 3 ) 

m=0 a=l 

To solve this equation it is convenient to remember that, as the Hilbert space 
TL n is a tensor product of a collective syndrome space S and a collective logical 
space A, we can always write the density matrix p as 

N 

p= £ \ m )LL{m'\ (g> Pmm>, (3.4) 
m,m'=0 

where the operators p mm ' act on the 2 fc ~dimensional space A. For p to be hermi- 
tian and normalized, the following identities must be satisfied: i) p mm i = Pm'm 
and ii) J2 m {pmm) = 1- Using thi s no tation one can rewrite both the corrective 
and the error-generating terms in (|3.3| ) as 

N N 

£ RmpRl = \0)ll(0\ ® £ U mPmm U m , (3.5) 

m=0 m=0 
N N N 

£ E aP E\ = £ £ |m + a) LL (m' + a\ ® e^-*^U aPmm iXJl, (3.6) 

a=l a=l mm'=0 

where errors are labeled using the commutation pattern c a (i.e. we choose to 
take a = c a ). These equations show that neither the erroneous nor the corrective 
quantum jumps couple diagonal operators (m = m') with off diagonal ones (m ^ 
m'). Using them, equations for diagonal operators p m = p mm can be easily found. 
To explicitly write down the equations for p m we will assume that the Hamiltonian 
Hq does not connect states with different syndromes. For this type of Hamiltonian 
(which are the ones required for a fault tolerant computation) we have Hq = 
J2m=o \m)(m\ (g) h m , where h m is a Hamiltonian acting in the logical subspace 
A. Using this, the master equation which includes the effects of Hamiltonian 
evolution, errors and error correction, reads: 

p m = -i [hm, p m ] ~ (t(1 ~ <W)) + N^')p m + 
N 

+ £ (tmOlUlpJJa + iU aPm+a U}) . (3.7) 
a=l 

As we can see from (|3.7|), all operators tend to be damped by the corrective term 
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except pq. On the other hand, errors tend to induce transitions between different 
syndromes at a rate proportional to 7'. From these equations we can first compute 
the probability for detecting the m-th syndrome, which is nothing but the trace 
of p m : defining p m = Tr(p m ) and taking the trace in we find: 

N 

Pm = -(7(1 - *mo) + N^')p m + ^ (drnOlPa + 1 Pm+a) • (3.8) 

0=1 

Solving this system, we obtain the probability for the correct syndrome 

»<" - 7 + (n + by + ("■"" - W(WT VP ) exp( - A "' (3 ' 9) 

where the decay rate is A = 7+ (N + l)^'. This rather simple formula, valid for all 
stabilizer codes in the limit of continuous correction, shows that the probability 
for detecting the correct syndrome tends to be of order unity if: i) the corrective 
coupling constant dominates over 7' (7 ^> 7' defines the strong correction limit) 
and ii) times are large as compared to t c = I/7. 

To solve the master equation finding the operators p m we will be more specific 
restricting ourselves to consider codes protecting k = 1 qubit. In this case the 
operator p m acts on a 2-dimensional Hilbert space and can be written as a linear 
combination of Pauli matrices: 

Pm = 2 (PmX + r m d^j . (3.10) 

where a is the vector of Pauli matrices and the polarization vector f m is the 
expectation value of a in the state p m (note that p m is not normalized since p m = 
Tr(p m ) is the probability for detecting the m-th syndrome). The master equation 
( |3.7| ) translates into a linear system of equations for the polarization vectors r m . 
Thus, the evolution of the quantum state is determined by the dynamics of the 
vectors f m , each one of which lives in its own Bloch sphere. To obtain equations 
for r m we assume that the Hamiltonian h m is proportional to a Pauli matrix 
(which induces a rotation of r m on the Bloch sphere). Thus, if h m = f2<7fc/2 we 
find 

N 

r m = Q,xr m - (7(1 - 5 m0 ) + N^')r m + ^A a (^5 m ar a + 7'r m+Q ) . (3.11) 

0=1 

Here, Q = ilk and the matrix elements of the 3x3 matrices A a are: 

(A«)y = \k 3 Tr{a,\J a f. (3.12) 

It is worth noting that, as for stabilizer codes the operators U a are Pauli matrices, 
the only possible values of the diagonal elements A at a are ±1. The equations for 
r m can be analytically solved in some simple cases. Here, we consider only the code 
protecting k = 1 qubit against phase errors affecting any one of the n = 3 carriers. 
In this case we can use the results given in the previous section to show that the 
matrices Ai and A2 are equal to the identity while A3 = diag(—l, —1, +1). Using 
this, we can obtain the equation for the polarization vector tq (for convenience 
we avoid writing the subscript and denote r = r~o). It turns out that r is 
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coupled through ( 3.11[ ) with two other vectors: fj = r\ + r*2 (the sum of the two 



polarization vectors corresponding to syndromes for which the recovery operator 
is U = I) and fz = fz (the polarization vector corresponding to the syndrome for 
which the recovery operator is U = Z). Thus, one can show that ( |3. 11 ) reduces 
to: 

f = x r — 37 r + (7 + 7') (rj + fz) , (3.13) 
r\ = n X rj - (7 + 3r/)rj + 7' (2r + 2f z + A 3 rj) , (3.14) 
= x - (7 + 37')rz + 7' (fj + A 3 r) . (3.15) 

where = (f2,0, 0). Details of the solution of this system will be analyzed else- 
where (Paz & Zurek, 1997). Here we will discuss the simplest case where an 
exact solution is possible. In fact, when O = (free evolution), the above system 
reduces to: 

f = -37V+ (7 + Y)i? (3.16) 
f = _( 7 + 3 7 ')^ + 3 7 y + 27'A 3 w (3.17) 

where the vector u is defined as u = XIaLi ^-afa = + ^z- Solving this system 
we can find the three components of the polarization vector f = (x,y,z). If the 
initial state belongs to the m = syndrome subspace (i.e., p(0) = |0)(0| ® Po(0)) 
we find: 

z(t) = z 0P0 (t) = z + (l - l±fy exp(-Ai)) , (3.18) 

y ^ = yo \_-x + ~ A +) exp ( ~ A -^ + ( A - " 37 ') exp ( ~ A+ ^ ' (3 - 19) 

where the decay rates are A± = + 7/2 =F y M7' 2 + 477' + 7 2 /4 and x(i) obeys 
an equation which is identical to that of y(t). 

Before discussing these results it is worth introducing a measure of the effec- 
tiveness of the error correcting process. For this we will use the state fidelity, 
which can be defined as 

F = Tr (p id (t) ]T R^pifyR^j , (3.20) 

where p(t) is the solution of the full master equation and is the ideal den- 
sity matrix obtained by unitarily evolving the initial state using the appropriate 
Hamiltonian. For the simplest case we are considering here (Hq = 0) the fidelity 
is determined by the scalar product between the initial state f(0) and the vector 

(f + v)(t). Thus, a simple calculation shows that F = \ (l + |r(0) (r + v) (t) 

which, using the above results, can be written as 

F = \ ((l + zl) + (4 + yl) x l _ x (A_ exp(-A+t) - A + exp(-A_t))^ 

(3.21) 

A first comment on the results is that there are some states for which the fidelity 
is always maximal. In fact, if the initial state is completely polarized along the 
z-axis, the code will always correct with ideal fidelity. The reason for this has 
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nothing to do with continuous error correction but is a simple consequence of 
the fact that for such states all the errors act trivially on the collective logical 
space A. Only states which have non vanishing components on the (x, y)-plane 
are nontrivially affected by the errors we are considering here. 

It is useful to analyze our results in various limiting situations. For the simplest 
case with no continuous correction (7 = 0), the fidelity is 

F = i ((l + zl) + (xg + vl) \ (3exp(-27't) - exp(-6 7 't))) . (3.22) 

This expression shows that, without continuous correction, the fidelity decays at 
a rat e fixed by 7' (the initial decay is quadratic because the fidelity defined in 
( |3.20| ) includes a final correction event at time t, which is sufficient to assure that 
for short times F ~ 1 — 67 /2 i 2 ). The effect of continuous error correction is clearly 
seen when comparing the above expression with the one arising in the limit of 
strong correction (7 3> 7') where, to first order in e = 7Y7, we can approximate 
A_ « 7(1 + 8e + 0(e 2 )) and A+ w 127e 2 ). In this case we find: 

z(t) =2 (l-3e + 3eexp(-7t)), (3.23) 

y(t) = y ((1 - 3e) exp(-12e 2 7 t) + 3eexp(- 7 i)) , (3.24) 

F = \[(l + zt) + (xl + yl) x 

x ((l + 12e 2 ) exp(-12e 2 7t) - 12e 2 exp(- 7 t)) ] (3.25) 

Two different time scales appear in the strong correction limit: O n the one h and, 
the short time scale t c = I/7 (after which the last terms in are 



exponentially suppressed) is associated with the corrective jumps. Over this time 

scale the normalized state p = | (i + fa/po^j lies on the surface of the Bloch 

sphere and is effectively protected from errors. On the other hand, on the much 
longer time £2 = t c /12e 2 there is an overall exponential damping of the (x, y)~ 
components of the Bloch vector which moves the normalized state p towards the 
interior of the Bloch sphere. On this time scale the error correcting code ceases 
to be efficient: fidelity decreases and purity of the quantum state is lost. 

Finally, we can briefly discuss the case where Hq is nontrivial (fi 7^ 0). For 
the three qubit code a fault tolerant rotation of the encoded state can be easily 
achieved by using the Hamiltonian Hq = QX3/2. Thus, as the encoding operator 
is such that CX%C^ = X3, applying Hq = QX3 during an appropriate time in- 
terval (to make an exact 7r-pulse) we obtain an encoded bit-flip. The reason for 
this is that while the operator X3 commutes with the generators Mj (and there- 
fore does not change the syndrome) it anti-commutes with the logical pointer 
L = Z1Z2Z3. Therefore, applying the operator X3 is tantamount to an encoded 
X operator. With this Hamiltonian, the polarization vector r evolves according 
to the system ( 3.13| - j05 ) which can be numerically solved. As the structure of 



the equation is similar to the one found in the free case one expects the results 
for the fidelity to be essentially the same we described above. The only relevant 
exception is that the initial state with Zq = 1 no longer produces ideal fidelity. In 
Figure 2 we show the time dependence fidelity for such initial state and display 
the (y, z) components of the polarization vector f. 
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Figure 2. Fidelity and non vanishing components of the polarization vector with and without 
continuous error correction. The initial state is polarized along the z axis and the Hamiltonian 
is Hq = a x /2 (the strength of the erroneous jumps is also set to unity). 
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